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Ordering and quantum fluctuations of orbital degrees of freedom are studied theoretically for 
LaV03 in spin-C-type antiferromagnetic state. The effective Hamiltonian for the orbital pseudospin 
shows strong one-dimensional anisotropy due to the negative interference among various exchange 
processes. This significantly enhances the instability toward lattice distortions for the realistic 
estimate of the Jahn- Teller coupling by first-principle LDA+[7 calculations, instead of favoring the 
orbital singlet formation. This explains well the experimental results on the anisotropic optical 
spectra as well as the proximity of the two transition temperatures for spin and orbital orderings. 

PACS numbers: 72.20.-i, 71.70.Gm, 75.30.-m, 71.27.+a 



Orbital degrees of freedom are playing key roles in 
magnetic and charge transport properties of transition 
metal oxides jjj. Especially it has been recognized that, 
in these strongly-correlated systems, spatial shapes of or- 
bitals can give rise to an anisotropic electronic state even 
in the three-dimensional (3D) perovskite structure [§,1). 
There the spin ordering (SO) and orbital ordering (00) 
are determined self-consistently 

Perovskite vanadium oxides, AVO3 (A is rare-earth el- 
ement), are typical t2 g electron systems which show this 
interplay between orbital and spin degrees of freedom 
|^|,|U^J|,^,|lO| . Both magnetic and orbital transition tem- 
peratures, Tn and T Q , respectively, change systematically 
according to the ionic radius of the A atom which 
controls the bandwidth through the tilting of VOg octa- 
hedra. For smaller ionic radii (smaller bandwidth) such 
as A=Y, T Q for 00 of G-type (3D staggered) is much 
higher than T N for SO of C-type (rod-type) As 
the ionic radius increases (the bandwidth increases), T 
decreases while Tn increases, and finally they cross be- 
tween A=Pr and Ce pi]] . In LaV03, the SO occurs at 
Tn — 143K first, and at a few degrees below Tn the 00 
takes place [[[2) . A remarkable aspect here is its proxim- 
ity of Tn and T Q , which is also observed for all the com- 
pounds with T N > T , i.e., CeV0 3 Jll| and Lai^Sr^VOs 
(x < 0.17) @. Therefore, in LaV0 3 , the magnetic cor- 
relation appears to develop primarily and to induce the 
orbital transition immediately once the SO sets in. 

Another interesting aspect of LaV03 is the large 
anisotropy in the electronic state, which has recently 
been explored by the optical spectra Q . Figure |l| shows 
the temperature dependence of the spectral weights, I c 
along the c direction and I ao within the ab plane, which is 
obtained from the data in ref. . Here we define the spec- 
tral weight as an integration of the optical conductivity 
up to the isosbestic (equal-absorption) point at 2.8eV, 



namely, 1^ = |^ J Q c cr u {oS)duj (fi = c or ab), where 
mo and n are the free electron mass and the density of V 
atoms, respectively. The most striking feature is the tem- 
perature dependence. I c grows rapidly below Tn while 
I a b is almost temperature independent. Therefore the 
temperature dependence is almost ID although the ratio 
Ic/Iab ~ 2 is not so large. 
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FIG. 1. Temperature dependence of the spectral weight 
in LaVOa. Filled and open circles represent the data along 
the c axis I c and within the ab plane I a b, respectively. The 
inset shows the optical conductivity. The solid curves de- 
note the data along the c axis for T = 10, 142, 293K from 
top to bottom, respectively. The data within the ab plane 
are almost temperature independent which are shown by the 
dashed curves. 

Compared with e g systems, the Jahn- Teller (JT) cou- 
pling in tig systems is expected to be weak, and the quan- 
tum fluctuation and/or the singlet formation of the or- 
bital degrees of freedom is a keen issue. Recently, Khali- 
ullin et al. claimed that orbital singlet correlation along 
the c axis is the driving force to realize the ferromag- 
netic spin exchange in the C-type antiferromagnetic (AF) 
phase in LaV03 fL4|| . On the other hand, the C-type 
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SO state with the G-type 00 has been obtained by the 
mean-field theory Jl5| and the first-principle calculation 
which are justified for rather weakly-correlated cases 
and do not take account of quantum fluctuations seri- 
ously. There, the AF interactions within the ab plane 
concomitant with a single occupation of xy orbitals play 
a key role to stabilize the SO and 00. Hence, these two 
pictures are quite different. Since there are competing in- 
teractions with the orbital quantum nature, such as JT 
coupling and 3D orbital exchange couplings, it is highly 
nontrivial to what extent the quantum fluctuations are 
important under the realistic situations, nevertheless the 
quantitative study on this issue has been missing. 

In this Letter, we study ordering and fluctuations of 
orbital degrees of freedom in the C-type AF phase in 
LaV03. An effective orbital model including the JT and 
the relativistic spin-orbit couplings is derived using the 
parameters obtained by first-principle LDA+[7 calcula- 
tions and the optical experiments. We show that the 
model exhibits a strong ID anisotropy which explains 
well the experimental results for the optical spectra and 
the proximity of T to Tn . It is concluded that the en- 
hanced JT instability due to the ID confinement domi- 
nates the orbital singlet formation in LaV03. 

Now we derive the effective orbital model. We start 
from the strong-coupling limit of the Hubbard model 
with three-fold orbital degeneracy for the t2 g orbitals 
. The system contains two d electrons at each V atom 
which form the high-spin 5=1 state due to the Hund's- 
rule coupling. To focus on the orbital sector in this 
spin-orbital coupled Hamiltonian, we assume the C-type 
SO. This is a reasonable approximation because the spin 
5 = 1 has less quantum nature compared to 5 = 1/2 and 
the C-type SO is obtained by the mean-field calculation 
JlU and the first-principle calculation (l(| as mentioned 
above. At the same time, we assume that the xy orbital 
is singly occupied at each V atom, which drives the AF 
coupling in the ab plane ||l5| , |l6| . Then the second elec- 
tron goes to either yz or zx orbital. We assign a pseu- 
dospin state t z — ±1/2 for the occupancy of these two 
orbitals |[4| . Finally, our total Hamiltonian is written as 
H = Ti-orb + H"fb + Wjt +7~Lls, each term of which is 
described below. 

The orbital exchange term for the c direction is 

«Srb=X>c(f<-7S--l/4), (1) 

m 

while that within the ab plane is 

<b = E ( ~ J aV4 - J%/4) ■ (2) 

m 

Here the summations are taken for nearest-neighbor 
pairs. The coupling constants are given by 
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J ab - 3(f// _ Jh) + 3(f// + 2Jh) ± U + 2 J li ± U' (4) 
J'2 = U 2 /(U + 2J Y{ )+U 2 /U, 



(5) 



where U, U' and Jh are the intra-orbital, the inter-orbital 
Coulomb interaction and the Hund's-rule coupling, re- 
spectively. Neglecting the small tilting of VOg octahedra 
in LaV03, the transfer integral t is taken to be diagonal 
which strongly depends on the direction, for instance, 
i« = t$T = t and otherwise zero in the c direction. From 
the analysis in ref. , the parameters are estimated as 
U ~ 2.25eV, U' ~ 1.93eV and J H ~ 0.16eV. The trans- 
fer integral t is set to be 0.1 2eV based on the estimate of 
the bandwidth ~ leV in first-principle calculations [fl6|| . 
Then the orbital exchange interaction along the c axis 
is estimated as J c ~ 33meV, while J~ b ~ 2meV. We set 
J c = 1 as an energy unit in the following calculations. 

Here we point out two important features in eqs. (|l|) 
and (||). One is that the exchange in the c direction is 
Heisenberg-type while that within the ab plane is Ising- 
type. From this, one might expect a strong quantum 
fluctuation in the c direction as pointed out in ref. . How- 
ever, this quantum nature becomes relevant only when 
the JT coupling is negligibly small, and this is not the 
case in LaV03 as discussed in the following. The other 
important feature is the large ID anisotropy in the orbital 
exchange couplings. Note that the negative interference 
among different perturbation processes occurs in the in- 
plane coupling J~ b in cq. (Q), which results in the ratio 
of the exchange couplings J c / J~ b ~ 17. 

The JT coupling in the subspace of r is given by 



(0) 



J c = At A /(U' - J H ) 



(3) 



where Qi is the JT phonon coordinate at site i. We ne- 
glect the kinetic energy of phonons, namely, regard Qi as 
a classical variable. It is crucial to estimate the coupling 
constant g, and we have done the following first-principle 
calculation fTij] . Assuming the tetragonal symmetry, we 
calculate the total energy as a function of the JT distor- 
tion as shown in Fig. ||. This gives the JT stabilization 
energy ~ 27meV, which approximately corresponds to 
£jt = g 2 /8 in our model. Thus we obtain the estimate 
Ejt ~ 0.8J C (g ~ 2.6), which is appreciable and cannot 
be neglected as in ref. . 

The last term is the relativistic spin-orbit coupling, 
which may be important in the ti g systems. We obtain 
the effective Hamiltonian by projecting the original form 
7~Cls — Yli M-'i ' &i to the subspace of r by using the ex- 
perimental fact that the spins lie within the ab plane |L8| . 
Here Li is the orbital angular momentum. Since L x has 
matrix elements between xy and {yz, zx) in this case, the 
spin-orbit interaction is represented by Hls = J2iC T i> 
where £ = X 2 / A and A is the energy separation between 
xy and (yz, zx) orbital levels. This indicates that the 
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spin-orbit coupling corresponds to the pseudo magnetic 
field along the z direction. Using A ~ 20meV in V atom 
and the estimate of A ~ leV in the band calculation Jl6| , 
we estimate ( ~ 0.4meV^ 0.01J C . This is small enough 
to be neglected in the following calculations. 



0.03 
0.02 
0.01 



^tot o.OO 
(eV) 

-0.01 
-0.02 





b a ( 

V V / 







-0.03 

0.000 0.005 0.010 0.015 0.020 

S/(v-v) cubic 

FIG. 2. Total energy per V atom as a function of the 
Jahn- Teller distortion obtained by the LDA+U calculation. 
The inset shows the tetragonal distortion of V06 octahedra 
considered here. The curve is a fit by a quadratic polynomial. 

First we discuss the anisotropic temperature depen- 
dence in the spectral weight. The spectral weight is gen- 
erally given by the kinetic energy of the system in the 
ground state K as I = —ttK/AN where N is the system 
size and we set e 2 = 1 jl9|. In the strong-coupling limit, 
this is calculated by the exchange energy as shown be- 
low. Let us consider the Heisenberg model for the strong- 
coupling limit of the Hubbard model at half filling. The 
kinetic energy in the Hubbard model is calculated as K = 
td(H)/dt = d{H)/d\nt. Here the bracket denotes the ex- 
pectation value in the ground state. Since the exchange 
coupling J in the Heisenberg model is proportional to t 2 , 
we obtain/ = -{it / AN)x2d(H) / d\n J = -7r(H Kcis ) /2N 
where 7iHois is the Heisenberg Hamiltonian. Using this 
formula, we calculate the spectral weight for the present 
model by the exact-diagonalization in the ground state 
for 4 x 4-site lattice embedded in the ac plane. In Fig. || 
(a), the obtained values of AI C = — 7tJ c ^(t; • fj)/2N 
and AI a b = —^J~ b ^2(T^Tj)/2N are plotted as a func- 
tion of the JT coupling g. Note that these exchange cor- 
relations correspond to the enhancements of the spectral 
weights from the high-temperature limit to the ground 
state. The results indicate that, the contribution in the 
c axis is much larger than that in the ab plane. The 
ratio becomes larger than 20 for the realistic value of 
g ~ 2.6. This explains well the anisotropic temperature 
dependence of the spectral weig ht in Fig. 0. 

We also discuss the magnitude of the total spectral 
weight in the ground state. Figure |3| (b) shows the calcu- 
lated spectral weights = — n(H^ lb )/2N which include 
the contributions from the constants in eqs. ([!]) and (||). 
As shown in Fig. || (b), the ratio I c /I a b is almost 1 for 
the realistic value of g, which is much smaller than that 
of the temperature dependent part in Fig. || (a). This 



ratio is comparable but smaller than the experimental 
result Ic/Iab ~ 2 in Fig. |l|. This might be due to the 
facts that (i) the cut-off energy 2.8eV for the integration 
is not large enough to take all the contribution of I a b, 
and (ii) the transfer integrals for the different orbitals 
are different because they depend on the energy differ- 
ence Ad p between the d and the oxygen p orbitals as the 
energy denominator. Since the huge anisotropy in tem- 
perature dependences is reproduced with the moderate 
anisotropy in the total weights, we believe that the or- 
bital ID confinement in our model plays a major role in 
the anisotropic electronic state in this material. 
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FIG. 3. (a) Contributions depending on temperature and 
(b) total weights in the ground state for the spectral weight 
as a function of the Jahn- Teller coupling g. 

We now turn to the discussion on the transition tem- 
peratures. In order to obtain the phase diagram at finite 
temperatures, we apply the transfer matrix method pQ ] 
combined with the unrestricted Hartree-Fock approxima- 
tion for the interchain coupling J~ b . We start from a 
configuration of {Qi} and calculate the expectation value 
(rf ) numerically. Then, from the self-consistent equation 
Qi = —g{rf) which is obtained as the energy minimiza- 
tion for Qi, we define a new configuration of {Qi}. We 
repeat this procedure until {Qi} is optimized and the en- 
ergy is minimized. In this treatment, J~ b is effectively ab- 
sorbed into the JT coupling g as ^ i {gQi — 2J~ b (T^))T^ = 
~ J2i 9 2 { T i) T i where g 2 = g 2 + 2J~ b . The obtained tran- 
sition temperature for the orbital/lattice ordering tran- 
sition is plotted as a function of g in Fig. |]. We note 
that for large g, the transition temperature diverges as 
T ~ g 2 /4 which corresponds to the JT energy gap. How- 
ever, this is an artifact of the mean-field-type treatment, 
and in reality T Q should stay at a constant of the order 
of J c since the coupling by J c between the neighboring 
sites determines T Q in the limit of large g. 

Considering the realistic value of J c ~ 33meV and 
Sjt ^ 27meV {g ~ g ~ 2.6) for LaV0 3 , the orbital tran- 
sition temperature is estimated as T Q ~ 800K which is 
much higher than the observed Tn — 143K as indicated 
in Fig. |4j One might think that this contradicts with 
the experimental fact Tn > T° xp . However the phase 
diagram is obtained assuming the C-type SO, which in- 
duces the ID confinement of the orbital degrees of free- 
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dom with the enhanced J c . Note that the disorder of the 
spins should reduce the effective orbital exchange as eas- 
ily shown in the spin-orbital coupled Hamiltonian. ft is 
also well-known that ID systems have an enhanced insta- 
bility to lattice distortions compared with higher dimen- 
sions. Therefore, T Q 1D = T Q under this ID orbital con- 
finement can be higher than T 3D without the SO when 
the JT coupling governs the 00 transitions. In real ma- 
terials with T a > In, T 3D decreases as the bandwidth 
increases, which indicates the relevance of the JT cou- 
pling. (If the 3D orbital exchange couplings dominate, 
T 3D should increase as Tn does.) Then, when the in- 
equality T 3D < T N < T D 1D is satisfied, the 00 transition 
with the JT lattice distortion should take place as soon 
as the SO grows and induces the ID confinement in the 
orbital channel. In this scenario, comparing with 
Tn, we can estimate the lower bound for the value of g 
to realize this proximity of the transition temperatures. 
Our estimation for this lower bound from Fig. [| is g ~ 1 
which is consistent with the estimate in Fig. ^[ 
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FIG. 4. Transition temperature of the orbital ordering un- 
der the ID confinement by the C-type AF spin order. For 
comparison, the experimental value of Tn in LaV03 is shown 
by the dashed line. The gray curve shows the expected Tq XP . 
See the text for details. 

Let us discuss this proximity of Tn and T° xp by 
the Ginzburg-Landau type argument. In this spin- 
orbital coupled system, we have the term in which the 
SO parameter M and the 00 one O are coupled as 
{a(T) - bM 2 }0 2 . Here a(T) is the coefficient for the 
second-order term for 00 without SO which is given 
by a(T) = a'(T - T Q 3D ). Our results indicate the re- 
lation a(T) - 6M s 2 at = a'(T - T G 1D ), where M sat is the 
saturated magnetic moment. T° xp is given by solv- 
ing the equation a(T |*P) - bM 2 (T ° xp ) = 0. Assum- 
ing M(T) = M satv /(Tv -T)/T N for simplicity, the 
difference between Tn and T° xp is given by ST/T^ = 
(T N - T cxp )/T N = (T N - T D 3D )/(T N + AT), where AT = 
T Q 1D — T 3D . Considering the systematic changes of Tn and 
T Q for A-site ions , we expect that T 3D is slightly lower 
than T N in LaV0 3 and CeV0 3 . Assuming T 3D = 0.8T N , 
we plot the expected T° xp in Fig. |^ as the gray curve. 
For the realistic value of g, we have T£ xp quite close to 



Tn as observed in these compounds. 

To summarize, we have investigated the role of or- 
bitals to understand the electronic state in LaV03. We 
have derived the effective orbital model with strong one- 
dimensional anisotropy assuming the C-type spin order- 
ing. We conclude that with the realistic Jahn- Teller cou- 
pling, the orbital ID confinement leads to an enhanced 
instability toward lattice distortions suppressing the or- 
bital quantum nature. This gives a comprehensive de- 
scription of the anisotropy in the optical spectra and the 
proximity of the critical temperatures of magnetic and 
orbital transitions. 
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